function grad = func(t,y)
tf = 1;
grad = zeros(2,1);
grad(1) = cos(atan(tf -t)) + y(2);
grad(2) = sin(atan(tf -t));


end









